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ERROR DIFFUSION WITH AVERAGED 
DIRECTIONAL BIASES 

CROSS-REFERENCE TO RELATED APPLICATION 

The subject matter of this appHcation is related to that of co-pending 
U.S. Patent Application Serial No. 09/593,445 filed June 15, 2000, by J. M. 
Gambaudo and C. Tresser for "Method and System Using Renormalized 
Pixels for Public Key and Compressed Images Watermarks on Prints" (IBM 
Docket YOR919990442US1) and assigned to a common assignee herewith. 
The disclosure of appHcation Serial No. 09/593,445 is incorporated herein by 
reference. 

DESCRIPTION 
BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention generally relates to digital printing and, more 
particularly, to digital printing methods using error diffusion (ED) but without 
the usual problems of anisotropy associated with ED methods. 

Background Description 

Error diffusion, as first described in "An Adaptive Algorithm for 
20 Spatial Greyscale" by R. W. Floyd and L. Steinberg in Proceeding of the SID 
17/2, (1976) pp. 75-77, is a popular technique for halftoning used today and is 
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considered to be one of the best, and the best among the techniques with 
similar operating time. This technique and others are reviewed in the book 
Digital Halftoning, MET Press, Cambridge, MA 1987, by R. Ulichney which 
is a general reference for digital halftoning. H. Kang, Color Technology for 
Electronic Imaging Devices, SPIE Optical Engineering Press, Bellingham 
Washington (1997) is another general reference on digital printing, including 
digital color printing. 

See also, J. F. Jarvis, C. N. Judice, W. H. Ninke, "A Survey of 
Techniques for the Display of Continuous Tone Pictures on Bilevel Displays", 
Computer Graphics and Image Processing, 5, pp. 13-40 (1976), P. Stucki, 
"MECCA-A Multiple Error Correction Computation Algorithm for BiLevel 
Image Hardcopy Reproduction", BM Res. Rep RZ1060 (1981), for 
introductory papers on Black and White (BW) ED that is designed for 
greysclale images. Color images can then treated using such ED on the 
fundamental colors Cyan (C), Magenta (M), Yellow (Y), and optionally Black 
(K), since K = C + M + Y, but using a separate K is more economical and 
usually provides a sharper black. 

U.S. Patent No. 5,070,413 to J. R. Sullivan, R. L. Miller, and T. J. 
Wetzel for "Color Digital Halftoning with Vector Error Diffusion" discloses 
ED working directly in the MX color space (one sometimes speaks of Vector 
Error Diffiision (VED) in this case). 

U.S. Patent No. 6,101,001 to C. P. Tresser and C. W. Wu for "Target 
Pattems Controlled Error Management" provides a description of a multitude 
of different algorithms of Error Diffiision. 

See J. P. Allebach and B. E. Rogowitz, eds.. Proceedings, SPIE - The 
International Society for Optical Engineering: Human Vision, Visual 
Processing, and Digital Display IV, vol. 1913, (San Jose, California), SPIE, 
Feb. 1993, for a collection of reprints that cover digital halftoning until 1993. 
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See also, U.S. Patent No. 6,614,556 to D.-E. Hong, C.-W. Kim, and 
G.-N. Boo for "Apparatus for Quantizing a Digital Image by Using an Error 
Diffusion Coefficient and Threshold Modulation in Zigzag Quantization", 
U.S. Patent No. 6,333,793 to K. Kobayashi for "Image Quality in Error 
5 Diffusion Scheme", U.S. Pat. No. 6,307,647 to A. C. Cheung, S. M. 

Heydinger, and S. T. Love for "Digital Halftoning With Error Diffusion", U.S. 
Patent No. 6,160,921 to G. G. Marcu for "Error Diffusion with Homogeneous 
Distribution in HighHght and Shadow Regions", U.S. Patent No. 5,684,932 to 
J. S. Shu for "Method and Apparatus for Dither Array Generation to Reduce 

10 Artifacts in Halftoned Image Data UtiUzing Ink Reduction Processing", U.S. 
Patent No. 5,668,638 to K. T. Knox for "Error Diffusion Method with 
Symmetric Enhancement", U.S. Patent No. 5,592,592 to J. S. Shu for "Method 
and Apparatus for Minimizing Artifacts in Images Produced by Error 
Diffusion Halftoning Utilizing Ink Reduction Processing", U.S. Patent No. 

15 5,521,989 to Z. Fan for "Balanced Error Diffusion System", U.S. Patent No. 
5,467,201 to Z. Fan for "Iterative Error Dififtision System",U.S. Patent No. 
5,051,844 to J. R. Sullivan for "Digital Halftoning with Error Diffusion", and 
U.S. Patent No. 4,654,721 to G. Goertzel and G. R. Thompson for "System for 
Reproducing Multi-Level Digital Images on a Bi-Level Printer of Fixed Dot 

20 Size", Qing Yu and Kevin J. Parker, "Stochastic Screen Design using 

Symmetric Error Compensation", University of Rochester (1997), and Ashish 
Jagmohan, Anshul Sehgal and Narendra Ahuja, "Isotropic Error Diffusion 
Halftoning", University of Illinois, for a collection of attempts at solving the 
problem of anisotropy in ED, with partial results. 

25 The problem of the consequences of anisotropy in ED (in particular 

the consequences on image quality) is considered as still mostly open despite 
the significant progresses which have been made since ED was invented, and 
as represented by the cited prior art. 
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Although considered a high quahty method, ED is marred with the 
problem that it is strongly anisotropic as the result depends strongly on the 
raster utilized to enumerate the pixels on a page. ED generates worm-like 
patterns, and other artifacts at the scale of the pixels. These artifacts affect the 
quality of the image, and are indeed attributed to the very anisotropic character 
of the algorithm. 

The importance of this invention is increased at a time when ink-jet 
printers arrive in the high end market and (the computational and data 
handling parts of the) processing times are bound to decrease dramatically for 
a few more years according to the perspective one has on the progress of core 
Information Technologies (IT), a decrease that one expects will be faster than 
the decrease in the mechanical part of the processing (actual printing and 
paper handling). 

SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide a cure for 
the anisotropy problem for multi-bit printers, that can print any of 2^ grey 
levels at each pixel for b bits, and more generally printers that can print any 
number G>1 of grey levels. 

According to the invention, there is provided a method to perform 
digital printing using error diffusion but without the usual problems of the 
anisotropy associated to that method. More particularly, the problems of the 
anisotropy associated to error diffusion are avoided, for printers that can print 
more grey levels than just white and black at each pixel, by choosing several 
rasters and error diffusions methods along these rasters (that may or not be all 
the same on the different rasters). Digital images are created for these pairs 
made by a raster and an algorithm followed by performing weighted averaging 



relatively low number of grey levels or colors at each pixel (for screen) or pel 
(for printing). We use the word pixel which is better known among non 
professionals as meaning either pel or pixel depending on the context. 

A multi-bit printer, that can print any of 2** grey levels at each pixel for 
5 b bits, and more generally printers that can print any number G>1 of grey 
levels, can be simulated by a one-bit printer (which can print only two grey 
levels, White or 0 and Black or 1, at each pixel) by creating super-pixels that 
are blocks of 2^- 1 (or more generally, B = G- 1) basic pixels. Figure la shows 
examples of two types of super-pixels with 2 bits (3 pixels per block) that are 
1 0 used in conjunction to pave the page. The representation of Va and of % can be 
fixed as for instance in Figure la, which might be most adequate for some 
printers (or other display systems) whose pixels are rectangular rather than 
being square. The representation of Vb and of % can instead be variable, using 
all or part of the choices available to represent Va and % as shown in Figure 
15 lb. Such choice of the representation when the grey level of the super-pixel 

has been determined to be Va or % can be made in many ways, for instance: 
• at random, 

or using a preselected ordering of representations to get the best 
possible image quahty: 
20 • out of a given printer (e.g., by making tests of various periodic 

sequences of the different renderings of a given grey value), 
or out of a family of printers (e.g., again by making tests of 
various periodic sequences of the different renderings of a 
given grey value), 
25 • the preselected ordering being made so as to either: 

depend on the grey level being represented (this is again done 
by laborious but easy to make testing), 
• or to not depend on the grey level being represented. 
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• or the choices of representation of Va and % to encode and hide some 
information in the super-pixels using the teaching of co-pending 
application Serial No. 09/593,445. 

Notice that similar procedures to what have been illustrated in Figures 
5 la and lb with two bits (or equivalently G = 4) could be used with other 

values of G, and one could expect small values of G to be more adequate for 
most displays. Also notice that the grey level values associated with groups of 
pixels have been computed in the example above by assuming that the printed 
pixels are perfect squares. For a lot of printers, there are in fact corrections to 

10 be taken into accoimt to get the exact effect of depositing some dots at some 
pixels. A detailed method to take account of the effective characteristics of 
dots is presented in U.S. Patent No. 5,946,455 to Tresser and Wu. It should be 
a simple matter for anyone well versed in the art of digital printing and finite 
mathematics to adapt such a method to the present invention. 

15 As well known in the art, once a problem of the kind considered in the 

present invention is solved for BW, the solution for color follows 
immediately. In fact, in the case of color printing, the invention can be used 
either by treating the fundamental colors independently, or by using VED. 
These considerations about the application for color to what has been 

20 developed for Black and White printing apply to both the cases of one bit 

printers (when super-pixels must be used to implement the present invention), 
or printers with G>2 grey levels, and now G>2 color intensities per color 
dimension. 

In fact, if one uses the CMY representation of the color space, where 
25 the ideal color space is the unit cube with axes corresponding to Cyan (C), 
Magenta (M), and Yellow (Y), White appear as the vector (0,0,0), and Black 
(K) as the vector (1,1,1). What is described along the grey axis in the Black 
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and White case is done along each of the coordinates and one sums up the 
results. One can also extract the black out of all pixels and treat that as a Black 
and White picture when a separate ink or toner is used for the Black for 
economy or better rendering of the Black; then the rest is treated along the C, 
M, and Y axes. One can also use VED as mentioned previously where the 
modified input is approximated by the closest comer of the color space and 
the error is treated as a three dimensional vector (when using VED, one can 
still separate the Black and White part). Instead of the CMY space, one can 
use more realistic color spaces where the simple cube geometry disappears 
which makes it hard to work one color direction at a time, but which allows as 
easy use of VED (Vector Error Diffusion) as the CMY space. 

Consider a printer that can print different G grey levels with G>2 at 
each pixel (not necessarily supposed to be uniformly distributed) on the grey 
levels scales (where 0 usually represents white and 1 represents black). We 
will choose some n with 0<n<G and then choose: 

n algorithms that may be all different, be all the same, or any choice in 

between, 

n different rasters, for instance, in the case when n=2, one could 
choose: 

• the standard raster, defined by going down the page line by 
line, and each line fi-om left to right, starting at the top left pixel 
which we designate by the pair formed by its column and line 
coordinates (1,1), and 

• the inverted standard raster, defined by going up the page 
line by line, and each line fi-om right to left, starting at the 
bottom right pixel (C,L) on a page with C columns and L lines 
of pixels. 
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Each ED algorithm Alg(M) in some library of ED algorithms can be 
defined by a set of weights defined for instance for the standard raster. As well 
know in the art, it is easy computation to adopt such weights to any raster, so 
that with: 

5 • L standing for the number of lines, 

C standing for the number of columns, 

R K, L. c standing for the K* raster choice, adapted to L lines and C 
columns, i.e., the pair 

Rk.l,c^(Ik.l,c» Fk.i,c) 

10 where: 

Ik,l,c = (h(0, K, L, C), v(0, K, L, C)), 
stands for the initial condition expressed by the original (starting) pixel of the 
raster, and where F^^^c describes the evolution of pixels positions according 
to that K"" raster as expressed by 
15 (h(m+l, K), v(m+l, K)) = ¥^^^1 (h(m, K), v(m. K))]. 

With now (a,b) standing for the vector fi-om the pixel being treated to the pixel 
whose error is being weighted, 

A(a,b, Alg(M),RK.i^c) 
will stand for the set of weights for the M* choice of an algorithm adapted to 
20 the K* choice of a raster. 

With such a set of weights A(a, b, Alg(i), Fj ^c). where now i is an 
index that run from 1 to n, (so that a pair (K, M) of a raster and an algorithm 
has been associated to i) and an input image Ijnp, it is well known in the art 
how to generate a corresponding digital output image I j j g (i) with a grey level 
25 at each pixel defined using M(G, n, i) grey levels using the ED precisely 
described by the chosen set of weights A(a, b, Alg(i), F; ^c)- Indeed, the 
modified input M j (h, v) at pixel (h, v) reads: 

Mi(h, V) = Ii„p(h, v)+ S(,,i,)i„v(i)W(.,b)Eri(h+a, v+b), 
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where V(i) is the set of vectors (defined by Alg(i) and Ri,L,c) where the errors 
needs to be taken into account to modify the input Ii„p(h, v) at (h, v), and 
W(a^b) is the weight at (h+a, v+b) in the set of weights A(a, b, Alg(i), F j ^c)- 
One can use errors zero out of the L lines and C columns to treat pixels near 
the boundary by the same formulas as the others. Other methods are easily 
devised; this is all part of the art of ED. 
Given Mi (h, v): 

the output I d i g (i)(h, v) at (h, v) for the ED algorithm Alg(i) is 

Clos(M(G,n,i)), i.e., the closest grey level g(i,h,v) among the M(G,n,i) 

ones that are available, and 

the new error Er j (h, v) computed at (h, v) is just: 
Er,(h, v) = M,(h, v)-I,jg(i)(h,v). 

The number M(G,n,i) is, as mentioned, the number of possible output 
grey levels (including 0 and 1), so that M(G, n, i)>l, and for the simphcity of 
the exposition we will assume that M(G, n, i) does not depend on i, so that we 
will use M(G, n) instead of M(G, n, i) to designate the number of grey levels 
that are possible at each pixel by each of the algorithms Alg(i). The number 
M(G, n) is not a (well defined) function as what counts is the number G of 
outputs after averaging according to the present invention as described below, 
which has to correspond to the printer being used, or to the number of grey 
levels that can be rendered at each pixel or super pixel. The values of the grey 
levels coming out from the algorithm determines the number of distinct 
outputs as soon as both n and G are large enough. For instance, when n=2 and 
the weights for the image averaging are Vi and 14, and for M = 4, the list of 
averaged outputs for outputs 0 = gl, g2, g3, g4 = 1 for both ED algorithms is: 
n g2 g3 
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e2 (g2^g3) g2+l 

g3 1^ 
^ ' 2 

1. 

One sees that according to whether g2+g3 = 1 or not, and according to 
5 whether ^^-^ = g3 or not, the number of averaged outputs that needs to be G 

is different. So for instance, for g2=y3 and g3=%, the Ust can be rewritten as: 
0, 1/6. V3, 14, 

10 1, 

with a list 0, i/e, Va, '/a, %, s/g,! of G = 7 different outputs, and we see that 
M(4, 2) = 4 is a possible choice. One can make M(G, n) a function if one 
imposes that the algorithms possible outputs are equally distributed among the 
possible M output grey levels and the weights for image average are all equal. 

1 5 Therefore, with averaging weights all set to be 1/n and with possible output 
grey levels for the ED algorithms all set to be of the form: 

0 _J 2_ (M-2) 

' (M-1)' (M-1)"" (M-1)' ' 

the number M(G, n) is a well defined function of G and n. Indeed, then the 
averaged outputs read 

20 0 — i (nM-n-1) 

' (n(M-l))' (n(M-l))' " ' (n(M-l)) ' ' 

which implies that G = n(M - 1)+1 = nM -n+1, so that M(G, n) is then given 
by 
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M(G. n) = i^lil^ 
n 

under all these simplifying assumptions (and only in such case). 

In the case of the so called multi-bit printers, the number G reads 2^, 
where b > 1 stands for the number of output bits (b=l would correspond to the 
5 single bit case, when only white or blacki can be rendered at each pixel), one has 
special forms 

M(2\ n) = 

n 

for the number of output grey level for the algorithms (Alg(i)). For instance, for 

b = 2,2*' = 4and^^li^ = l,for n = 2, 
n 2 

1 0 which cannot be as M must be an integer, 

(2^-l■^n) 6 . . 

-^^ = — = 2, for n = 3, 

n 3 

so that n = 3 is a possible choice, 

(2^-l^n) 7 p , 
= -,for n = 4, 

n 4 

which cannot be as M must be an integer, and no further value of n is possible, 
so that for b = 2 (or G = 4), n = 3 and M = 2 is the only possible choice under all 
the simplifying assumptions that we have made. The averaged outputs are then 

0,y3,%, 1, 

and the possible grey for Alg(l), Alg(2) and Alg(3) are just 0 and L 

According to the invention, the final image is obtained by averaging 
over the n images Idi g(i), but the raster may be chosen so that averaging can 
begin before the images I^i g(i) are computed so that less information needs to 
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be traded (another possible tradeoff between efficiency and quality). The 
averaging can be done with equal weights for all images, but one may prefer to 
use different weights for the different rasters. For instance, one can take two 
rasters, each providing an image with a two grey levels (0 and 1), and take one 
5 image for one of the rasters with weight %, and the other one with weight Va 
to get 4 grey levels (0, Va, %, and 1) for the averaged image. Altogether, there 
will need to be some tradeoff between the precision per image of the form 
Idi g(i) and the willingness one has of destroying the wormy artifacts andy any 
further anisotropy artifact of ED. The right number n will in general depend 

10 not only on the printer being utilized, but also on the image, with images being 
possibly classifiable according for instance to the violence of the contrast, the 
presence or not of lots of straight lines, etc., which may help choosing 
appropriate combinations of rasters, algorithms and weights. 

With reference now to Figure 2, one can follow how the invention can 

15 be implemented in a preferred embodiment. Block 5 is a library of ED 
algorithms, which can be kept for instance in standard form, i.e., by the 
weights that describe the algorithm with the standard raster which consists 
in following line by line, going down the line, and following each line from 
left to right. Block 10 is a library of rasters and of error diffusion algorithms 

20 (any given raster can come with any choice of error diffusion algorithm), 
where the raster can be for instance represented in the form 

where the index K designates the K* choice of raster, I l, c stands for an 
initial condition, i.e., the choice of an initial pixel on a page with L lines and C 
25 columns, F l, ^ stands for the function or lookup table that describes in which 

order the pixels are visited on the page to perform the chosen ED algorithm, 
so that with (h(0, K, L, C), v((0, K, L, C)) standing for the initial pixel of the 
raster R k. u c (for instance for the standard raster, say R j l, ^ » 
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(h(0, 1, L, C), v((0, 1, L, C)) = (1, 1), i.e., the top left pixel of the page for any 
L and C), and (h(m, K, L, C), v((m, K, L, C)) standing for the m* pixel along 
the K*^ raster, we have 

(h(m+l, K, L, C), v(m+l, K, L, C)) = F^^Mm, 1, L, C), v((m, 1, L, C)). 
For instance, with the standard raster 

(h(m+l, 1, L, C), v(m+l, 1, L, C)) = (h(m, 1, L, C)+l, v(m, 1, L, C)), 
ifh(m, 1,L, C)<C, 

(h(m+l, 1, L, C), v(m+l, 1, L, C)) = (1, v(m, 1, L, C)+l), 
if h(m+l, 1, L, C) = C, and v(m, 1, L, C) < L, and is not defined if 
h(m+l, 1, L, C) = C, and v(m, 1, L, C) = L, i.e., the raster stops there. 

At block 10 there is also a library of ED algorithms indexed by M. If 
the algorithm Alg(M) is chosen with the raster Rk, u c > its set of weights is 
adapted as the set 

{W(3,b)}=A(a,b,Alg(M),R^^e) 
such that the modified input from an input image Ij at pixel (h, v) reads: 

MK,M(h, V) = Ii„p(h,v)+ 2(a.b)inv(M) W(,,b)ErK,M(h+a, v+b), 
where V(M) is the neighborhood associated to the algorithm Alg(M) and the 
raster R l, c » and Er (i, j) stands for the error at pixel (i, j) for the choices 
K for the raster and M for the ED algorithm and whichever image is being 
treated. 

At database 20, there are sets of weights for averaging the output 
digital images Iji g(i). A set of weights may be chosen once and for all, or may 
depend on the rasters being used, the algorithm being used, and the input 
images of some classes of input images (such as dark versus clear, landscapes 
or portraits or neither, etc.). 

At input block 50 is the input image that is a matrix (or array) 
Iinp(h, v). We assume that G > 2 stands for the number of grey levels that can 
be rendered at each pixel (in the case when G = 2, one builds super-pixels as 



4 



15 

explained previously to have an effective number of grey level greater than 
two). At function block 100, one chooses some number n, with 0 < n < G, of 
rasters that will be used for the process of digital printing according to the 
present invention. This number n of rasters and of algorithms (although all 
5 algorithms may correspond to the same one on the standard raster) will be 

chosen at function block 120 from the library 10, i.e., for each i that satisfies 
the inequality 

(hmits on i) 0 < i < n+1 
a raster among the rasters R k, u c that we designate by R c »and 
10 • an algorithm among the algorithms Alg(M) that we designate by 
Alg(i). 

With i initialized to be 0 at function block 100, i is incremented at 
function block 110, replacing the input i by an output new value i+1. At 
function block 120, the raster Rj^c and the algorithm Alg(i) with order 

1 5 number i for the pass of the method, as described previously, are picked. 

There s a set of ED weights A(a, b, Alg(i), Ri,L.c) associated to the choices of 
the raster RjL,c and the algorithm Alg(i). At function block 120 is also chosen 
the number M(G, n) of grey levels among which outputs will be sought by the 
algorithms Alg(i), and the grey levels themselves among which the outputs 

20 will be sought. In fact, M(G, n) and the grey levels themselves among which 
the outputs will be sought are part of the specifications of the algorithms, 
besides the ED weights A(a, b, Alg(i), R^^^cl The number M=M(G, n) and 
the grey levels themselves among which the outputs will be sought are such 
that after averaging the possible outputs at each pixel, one gets the G grey 

25 level that are printable at each pixel by the printer (but the grey levels 

themselves among which the outputs will be sought are abstract in the sense 
that they only participate to average and need not be printable). So if the grey 
levels themselves among which the outputs will be sought are 
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0 = gl < g2 <...< gM = 1, (where we have ignored the freedom of dependence 
on i for an easier exposition, but that freedom can be used in the invention as 
anyone versed in the art of combinatorics would understand readily), and the 
weights for averaging from database 20 are W1,W2, ...,Wn, the possible 
outputs are 
0, 

Wl.g2,W2.g2,...,Wn.g2, 
Wl.g3,W2.g3, ...,Wn.g3, 



10 Wl.g(M-l), W2. g(M-l), Wn.g(M-l), 

Wl, W2, Wn, 

Wl.g2+W2.g2, Wl.g2+W3.g2, Wl.g2+Wn.g2, 
Wl.g2+W2.g3, Wl.g2+W3.g3, Wl.g2+Wn.g3, 



15 Wl.g2+W2, Wl.g2+W3, Wl.g2+Wn, 

Wl.g3+W2.g2, Wl.g3+W3.g2, Wl.g3+Wn.g2, 



and so on. For instance, in the case of two rasters with equal weights '/z, if the 

grey levels among which the outputs will be sought for individual images are 

20 0 = gl <g2 < ... < gM = 1 , the possible averages are 

Q g2 g3 g4 
' 2 ' 2 ' 2 

o2 (g2^g3) (g2^g4) (g2H-gl) 

„3 (g3-^g4) (g3+gS) (g3-Hgl) 
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g(M-2).(g(M-2)+g(M-l)). (g(M- 2)^-1) 

2 

g(M-l), (g(M-^^)-l) 

i: 

When M = 4, the list becomes 

0, g3 
2 2 



g2 (g2^g3) (g2-^gl) 
2 ' 2 



gs.ii^lil) 



1, 

so that M(7, 2) = 4 works with g2 = Va and g3 = % for seven outputs 0, i/e, Va, 
1 0 '/2, %, sk, 1 . When M = 3, the Ust becomes 

g2, 

1, 

so that M(5, 2) = 3 works with g2 = 'A for five outputs 0, 'A, 'A, When 
15 M = 2, the list becomes 

0,>/2 

1, 

which forces G = 3. 

As we have seen, and as well known in the art, the adaptation of an ED 
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to color once it is well defined for black and white poses no problem, for 
instance by using three independent color planes, or by first extracting the 
black part and then treating what remains with three independent color planes. 
At function block 130 the digital image Idig(i) corresponding to the 
5 input image I^^p at 50 and to the choices of Rj and Alg(i) and of the 
associated ED weights A(a, b, Alg(i), Ri,^c) at function block 120 is 
computed. This digital image Idig(i) is then sent to a temporary storage 140. 
From temporary storage 140, a comparison is made at decision block 150 to 
determine whether or not i = n. If not, i < n and one follows the path 1 80 to get 

10 back to the adder for incrementing i at 1 10, from where one continues as 

before. If yes, i = n and one follows the path 200 that prompts averaging of all 
the digital images Idig(i), i = 1, 2,..., n corresponding to the input images Ij^p 
at input block 50, and the choices at function block 120, and accumulated in 
the temporary storage at 140, using the weights prepared according to all the 

1 5 choices at database 20, resulting to the final output image I ^ ^ „ according to 
the present invention at 300. 

The averaging process can be done on parts of images rather than on 
fiill images to save some space, or to judge the choice made of the rasters and 
algorithm before the whole processing is done. 

20 With reference now to Figure 3, there is shown the typical situation 

where a paper image 610 is scanned at scanner 620 to generate a digital image, 
although the digital image can be the form in which the image to be printed is 
given. After processing of the image in the computer 630, the digital image 
associated to the input image according to the present invention is printed on 

25 the digital printer at 640. The method can be installed in the computer using 

different means such as a diskette, or using transfer through a network, 
including a global network such as the Internet, and updates such as new 
rasters or new algorithms can be installed in the same way. 
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While the invention has been described in terms of a single preferred 
embodiment, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the appended 
claims. 
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